The Eliashberg function of amorphous metals 
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A connection is proposed between the anomalous thermal transport properties of amorphous solids 
and the low-frequency behavior of the Eliashberg function. By means of a model calculation we 
show that the size and frequency dependence of the phonon mean-free-path that has been extracted 
from measurements of the thermal conductivity in amorphous solids leads to a sizeable linear region 
in the Eliashberg function at small frequencies. Quantitative comparison with recent experiments 
gives very good agreement. 
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The functional form of the Eliashberg function, 
a 2 F(w), in amorphous metals has been the subject of 
debate for a long time. Experimentally, there is uni- 
form agreement that in amorphous simple metals the 
low-frequency part of the Eliashberg function is strongly 
enhanced over what is observed in crystalline materials.!!! 
The functional form at low frequencies, although dif- 
ficult to determine experimentally, is usually found to 
be linear .□ This experimental observation has given rise 
to a substantial theoretical debate. BergmannH and 
many othersH argued that a disorder enhancement of the 
electron-phonon coupling leads to a linear low-frequency 
behavior of the Eliashberg function, in accord with ex- 
periment. Tiiis was disputed by SchmidEl and by Keck 
and SchmidQ, who claimed that this disorder enhance- 
ment was spurious, and due to an incorrect application 
of the Frohlich model to disordered materials. Schmid's 



calculation gave instead a 2 F(u 



0) 



in disagree- 



ment with_experiment. The cu 3 result was also confirmed 
by others .0 Since different models and different calcula- 
tional methods had been used by these various authors, 
the theoretical problem was widely considered unsolved 
for a long time. This disagreement was finally settled 
by Reizer and Sergeyev,H who explicitly pointed out the 
errors in Ref. |] and showed that a correct calculation 
leads to Schmid's result independent of the model and 
the method used. This led to the unsatisfying situa- 
tion that a theoretically credible result, viz. Schmid's 
a 2 F(uj — > 0) ~ u 3 , disagreed with experiment, while in- 
correct arguments led to a 2 F(uj — > 0) ~ u> in agreement 
with the experimental observations. This situation was 
summarized by one of usp who also argued that no exist- 
ing comparison between theory and experiment had been 
careful and accurate enough to really rule out Schmid's 
result. However, recent experiments by Watson and Nau- 
gle have shown that Schmid's result is not compatible 
with experiments on amorphous CuSn alloys.£3 

At this point it is important to distinguish between 
Schmid's general theory for the electron-phonon coupling 
in impure metals, and his and others' specific model cal- 
culations. Since the work of Reizer and Sergeyev showed 
that Schmid's general theory is physically correct, the 



search for reasons behind the discrepancy between his 
results and experiment should then turn to the model 
assumptions. On the electronic side the main assump- 
tion is that of nearly free electrons. It is hard to see how 
this could qualitatively fail in simple metals as long as the 
resistivities are moderate, and the electronic states are ef- 
fectively three-dimensional. On the phonon side, Schmid 
assumed undamped Debye phonons, an assumption that 
is necessary in order to obtain the cj 3 -law as the asymp- 
totic low-frequency behavior. As was shown in Ref. 
the inclusion of phonon damping by electrons leads to 
a linear low-frequency asymptotic behavior (albeit with 
a prefactor that is too small by several orders of mag- 
nitude to explain the experimental results), which then 
crosses over to Schmid's u 3 -law. The prefactor of the 
linear term is proportional to the phonon damping. This 
raises the possibility that very strong phonon damping 
(which would have to be of other than electronic origin) 
might lead to a linear term in a 2 F{uj) whose prefactor 
is large enough to account for the experimental observa- 
tions. 

In order to pursue this last point, let us recall that 
besides the problems with the Eliashberg function men- 
tioned above, amorphous materials have properties of en- 
tirely phononic origin that are hard to understand. In 
particular the thermal conductivity, K, shows an enig- 
matic behavior. Even though it has been stressed that 
the thermal conductivity is not understood in any tem- 
perature region, the general ph^juomenology is clear, con- 
sistent, and well documentedilM 2 ! As a function of tem- 
perature T, the thermal conductivity behaves like n ~ T 2 
for T/Q < 10~ 2 , with 9 the Debye temperature. The ori- 
gin for the phonon scattering in this region is not known 
for certain. The phenomenological two-lev&Lsystem con- 
cept has often been invoked in this context ,113 but no con- 
sensus has ever been reached. For 10~ 2 < T/Q < 1CP 1 
the thermal conductivity is approximately independent 
of T. This is the so-called plateau region, which is 
characterized by strong, and strongly frequency depen- 
dent, phonon scattering of uncertain origin. Finally, 
for T/Q > 10 _1 the thermal conductivity becomes In- 
dependent again, but it is not even clear whether the 
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heat transport in this region is by phonons, much less 
what the scattering mechanisms are. 

This poor state of physical understanding notwith- 
standing, the above phenomenology is remarkably uni- 
versal, and seems to bp-characteristic oLamorphous mate- 
rials, both insulatingji 1 ! and metallic.^ 3 ] It has been used 
to deduce the following behavior of the phonon mean- 
free path, l p h, as a function of frequency. For frequencies 
lu < 10^ 2 kB&/h, lph is a linear function of frequency. 
For intermediate frequencies, 10~ 2 < hco/ksQ ^ 10 _1 , 
lph goes as a high power, n, of frequency, n has been 
reported to be at least 4, and possibly larger. This in- 
termediate frequency regime corresponds to the plateau 
region in the thermal conductivity. At still higher fre- 
quencies, uj £ 10 _1 fcsG/?i, the phonon-^nean free path 
either becomes frequency independent Jlil or is a linear 
function of frequency agaiii.til 

In this paper we propose a connection between the 
thermal properties of amorphous materials as described 
above, and the low-frequency behavior of the Eliashberg 
function. In particular we show that Anderson's phe- 
nomenological functional form of the phonon mean-free 
path, if used in Schmid's theory for the electron- phonon 
coupling, explains the observed behavior of the Eliash- 
berg function as well as the observed behavior of the ther- 
mal transport. Let us start from the expression for the. 
Eliashberg function, based on Schmid's general theory,0 
that was derived in Ref. ||, 

(la) 

Here D b H (q, lu) is the retarded phonon propagator, whose 
imaginary part reads, 



Im£>£(q,w) 



4w^(q) 76 (q) 



lu^+Alu 2 ^) 



(lb) 



where 7(q) is the phonon damping coefficient. In writing 
Eqs. (1) we have assumed a free electron model with Np 
the electronic density of states per spin at the Fermi level. 
We have also assumed Debye phonons with one longitudi- 
nal and two transverse branches labeled by b (b = L,T), 
speed of sound q,, and dispersion w&(q) = c b q. a&(q) in 
Eq. (|la|) is the electronic contribution to the sound atten- 
uation-epefficient, for which we use the standard Pippard 
result Jia 



a b (q) = Kbfb(ql) 



(2a) 



where n b — {vp / Cb){p e / Pion) /I with vf the Fermi ve- 
locity, I the electronic mean-free path, and p e and pion 
the electronic and ionic mass density, respectively. The 
functions /t t are given byjia 



h(x) 



1 x 2 arctan(x) 
3 x — arctan(x) 



(2b) 



Jt{x) = — — 7 |"22: 3 + 3a;-3(a; 2 + l)arctan(a;)l . (2c) 
2x A 

With phonon damping exclusively by electrons, as was 
assumed in Ref. |9[ one has 7b(q) = Q,ab(q)/2. Here, 
however, we will consider the possibility of nonelectronic 
contributions to 7b (q). Accordingly, we write 



lb (q) =Jc b q D g(q/q l 



(3) 



where qo is the Debye wavenumber, 7 is a number, and g 
is some function that determines the wavenumber or fre- 
quency dependence of the phonon damping. The latter 
we model after Anderson's proposal,! 14 ] which has been ex- 
tracted phcnomenologically from thermal transport mea- 
surements in amorphous materials. Anderson's model 
consists of the three distinct regions mentioned above: 
(1) A low frequency region where the damping is a linear 
function of frequency, (2) an intermediate region where 
the damping goes as a large power of the frequency, and 
(3) a high frequency ifegion where the damping is inde- 
pendent of frequency.tS The intermediate region corre- 
sponds to the characteristic plateau that is observed in 
the T-dependent thermal conductivity. We thus model 
the function g(x) in Eq. (||) as, 



g(x) = 10 n y 



x/y + (x/y) n 
10™ + {x/y) n 



(4) 



Here y is the onset of the plateau region in units of the 
Debye wavenumber, the width of the plateau region has 
been assumed to be one decade, and n is the power that 
characterizes the frequency dependence of the phonon 
mean-free path in the plateau region. 

Before we turn to a numerical evaluation of the inte- 
gral, Eq. (|la|), that determines a 2 F(uj), let us consider 
the low-frequency behavior analytically. Asymptotically, 
a 2 F(uj) ~ uj/uj a , with a slope u>~ . The latter we esti- 
mate for a clean system, i.e. in the limit I — > 00. In this 
limit only longitudinal phonons contribute, and we can 
use the asymptotic form of the function /j, in Eq. (pb]), 
fh{x — * 00) = nx/6. Then we obtain, 



LU a 6lT k F \C L J Pion J a x 



(5) 



Typical parameter values are qnl^F ~ 1> vp/c^ « 10 3 , 



and p e l p. 
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For the parameters y and n in Eq. 



(|4|) we takeu-H y « 0.02 and n « 4. Finally, 7 determines 
the overall scale for the phonon mean- free path l p h- A 
typical value is l p h ~ 1cm at a frequency of 1GHz. With 
cl ~ 2 x 10 5 cm/s this corresponds to 7 ~ 2 x 10 -4 . This 
yields ep/uj a sa 1, 700. With a Fermi energy tp rj lOeV 
we obtain uj a « 6meV. This value for uj a is of the same 
order of magnitude asJ-he one typically obtained from 
tunneling experiments. ut3 

Now that we have seen that we obtain promising re- 
sults for a 2 F(ui — > 0) with reasonable parameter values, 
let us calculate a 2 F(ui) numerically, and comparp-quan- 
titatively with experiments. Watson and Naugldi3 have 



2 



performed a detailed study of amorphous SnCu. For the 
stoichiometry Sn.s7Cu.13 they quote the following param- 
eter values: ep = 1-54 x 10 _11 erg, kp = 1-59 x 10 8 cm _1 , 
v F = 1.84 x 10 8 cm/s, I = 9.58 x l(T 9 cm, q D = 1.31 x 
10 8 cm-\ c L = 1.6 x 10 5 cm/s, c T = 8.1 x 10 4 cm/s. Of 
these, the electronic parameters are much better known 
than the two sound velocities. Using these parameters, 
as well as n = 5, y = 0.015, and 7 = 8.0 x 10 -5 , we have 
calculated a 2 F(uj) for frequencies up to 1.6meV, which 
was the lower frequency cutoff in the experiment of Ref. 
|l0| . The high frequency behavior resulting from our cal- 
culation would not be realistic anyway due to our using a 
Debye model. The result was shown in Ref. |l0|, and was 
used as low-frequency input in a McMillan-Rowcll inver- 
sion procedure to obtain a 2 F from tunneling data. It is 
also shown again as the curve labeled n = 5 in Fig. [I]. For 
the inversion procedure an overall factor multiplying the 
calculated a 2 F was used as a fit parameter. The need for 
such an overall scale factor is not surprising, given our 
free electron model. The factor used for the best fit is 
equivalent to a deviation of the density of states in Eq. 
( p"a| ) from its free electron value by 14%. A comparison 
between the calculated and the measured tunneling den- 
sity of states then provides a measure of how well the 
low-frequency input describes the actual system. Wat- 
son and Naugle found that our calculated a 2 F does very 
well, although not quite as well as if one assumes a strictly 
linear low-frequency behavior. It should be stressed that 
our calculation used strictly the parameters as provided 
by the experimentalists, some of which are not known 
very accurately. Since the inversion procedure is quite 
involved no attempt was made to fine tune the parame- 
ters. 

This result shows that Schmid's theory with a phonon 
damping that accounts for the thermal transport proper- 
ties characteristic of amorphous metals gives good agree- 
ment between the calculated Eliashberg function and 
tunneling data. In contrast, the same theory with phonon 
damping by electrons onlv is not capable of explaining 
the experimental results. 

In addition to this comparison between theory and ex- 
periment, let us demonstrate the effects of some parame- 
ter changes on a 2 F. We consider the four results for o?F 
shown in Fig. [I]. The curve labeled n = 5 was obtained 
with the parameters as given above. The slight bulge 
in this curve results from the leveling off of the phonon 
mean-free path at the high-frequency end of the plateau 
region. This moderates the rapid increase of a 2 F at lower 
frequencies, which is due to the strong frequency depen- 
dence of the phonon mean-free path. With a weaker 
frequency dependence of the phonon mean-free path in 
the plateau region, i.e. a smaller exponent n in Eq. (|J), 
the initial slope of a 2 F is much smaller, and over the 
frequency range considered a 2 F shows a purely positive 
curvature. Conversely, a still larger exponent n leads to a 
purely negative curvature of a 2 F. This is demonstrated 
111 Fig. g. The curves with stronger curvature all led to 
substantially less good agreement with experiment than 



the one for n = 5. We have also considered the sensi- 
tivity of the result to the ratio of the longitudinal and 
transverse speeds of sound, which is not known very ac- 
curately. We have found only a very weak dependence 
of the functional form of a 2 F on this ratio in the region 
1.8 < cl/ct < 2.5. Finally, we have changed the damp- 
ing parameter, 7, with all other parameters held fixed. 
This was found to have a very similar effect to changing 
n, with a 2 F changing from negative to positive curva- 
ture as 7 is increased, or the phonon mean- free path at a 
reference frequency is decreased. The effect of changing 
7 by a factor of ten was roughly equivalent to changing n 
by one. For instance, with n = 4 and 7 = 8 x 10 -4 we ob- 
tained a curve that was hardly distinguishable from the 
one for n = 5 shown in Fig. p]. Generally, we found that 
with reasonable parameters for simple metals we need 
4 ^ n < 6 in order for our explanation of the behavior of 
a 2 F to be viable. 

In conclusion, we have shown that Schmid's theory of 
electron-phonon coupling in impure metals can account 
for the observed low- frequency behavior of the Eliashberg 
function in amorphous simple metals if one assumes a 
strong phonon damping consistent with the one extracted 
from measurements of the thermal conductivity. While 
the physics underlying the strong damping is not known, 
this observation unifies two seemingly unconnected, and 
separately mysterious, properties of amorphous materi- 
als. It suggests that strong phonon scattering is a very 
fundamental feature of the amorphous state, and that un- 
derstanding its origin would explain many different prop- 
erties of amorphous materials at once. 
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sions. This work was supported by the NSF under grant 
numbers DMR-92-09879 and DMR-90-19525. 
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FIG. 1. Results for the Eliashberg function a 2 F as a func- 
tion of energy or frequency for different values of the exponent 
n in Eq. (Q). All other parameters were held fixed at the val- 
ues given in the text. a 2 F has been normalized by its value 
at u> = 1.6meV, which is 0.436, 0.279, 0.208, and 0.200 for 
n — 6, 5, 4, and 3, respectively. 
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